Novel massless phase of Haldane-gap antiferromagnets in magnetic field 
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The behavior of Haldane-gap antiferromagnets in strong magnetic field is not universal. While 
the low-energy physics of the conventional fD spin-f Heisenberg model in its magnetized regime 
is described by one incommensurate soft mode, other systems with somewhat perturbed coupling 
constants can possess two characteristic soft modes in a certain range of the field strength. Such a 
too-component Lutinger liquid phase is realised above the massive Haldane-gap phase, and in general 
above any massive nonmagnetic phase, when the ground state exhibits short range incommensurate 
fluctuations already in the absence of the field. 
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Quantum mechanical systems possessing a spectral gap 
in their ground state are usually very robust: the gap can 
persist even if relatively large perturbations are added to 
the Hamiltonian. This is exactly the case for the one- 
dimensional spin-1 antiferromagnetic chain where the ex- 
istence of an energy gap, the Haldane gap, El is well doc- 
umented and understood for the Heisenberg model. In 
the more general SU(2) symmetric bilinear-biquadratic 
spin-1 model 
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in which (3 = yields the conventional Heisenberg model, 
the Haldane gap survives in the whole range — 1 < (3 < 1, 
thus there is a large region in the /3-space where the bi- 
quadratic term has seemingly no effect on the low-energy 
physics .□ The Haldane gap disappears in the integrable 
critical points [3 = ±1 beyond which new phases appear. 
For (3 < — 1 translation invariance is spontaneously bro- 
ken and dimerization occurs. In the region [3 > 1 the 
bilinear-biquadratic model is believed to remain gapless 
with soft modes at momenta k = 0,±2-7r/3. Since the 
system is one dimensional no long-range order exists in 
this phase either: the correlation functions decay alge- 
braically. 

The emergence of the Haldane gap in the vicinity of the 
f3 = ]_point, the so called Uimin-Lai-Sutherland (ULS1 
point J3 has been investigated recently by Itoi and Katocl 
in the absence of magnetic fields h = 0. The ULS model 
has an SU(3) symmetry which is broken down to SU(2) 
when [3 is tuned away from 1. They identified the criti- 
cal theory of the ULS point with the k = 1 SU(3) Wess- 
Zumino-Witten-Novikov (WZWN) model and concluded 
that the SU(3) symmetry breaking perturbation, repre- 
sented by the deviation term with — 1) in the Hamilto- 
nian, is irrelevant for (3 > 1, i.e., the system remains crit- 
ical there, but it becomes marginally relevant and gives 
rise to a dynamic mass generation (the Haldane gap) for 
(3 < 1. The transition is of the Berezinskii-Kosterlitz- 
Thoulcss (BKT) type, i.e., the gap open exponentially 



slowly away from the ULS point. Tiese findings are sup- 
ported by earlier numerical results. El 

When there is a magnetic field present the SU(2) sym- 
metry of the bilinear-biquadratic model breaks down to 
U(l). At the ULS point, however, where the symmetry 
is higher, the quantities N + ,Nq, and 7V_, denoting the 
numbers of +, 0, and — spin states in the wave function, 
with N + + Nq + N- = N the length of the chain, are in- 
dependently conserved. Thus the remaining continuous 
symmetry for (3 = 1, h > is U(l) x U(l). For (3 < 1 
when the strength of the field is higher then the value of 
the gap, the Haldane gap collapses and the low-energy 
physics of the magnetized system is governed by gapless 
excitations. The emerging periodicity is a function of the 
field strength and in the generic situation incommensu- 
rate. Finally, when the field is strong enough all the spins 
align and the magnetization saturates. 

At the ULS point, where the Haldane gap no longer ex=. 
ists, the analysis of the Bethe Ansatz equations showedcl 
that the magnetization growth is not smooth. There is 
a second order phase transition which leads to a cusp 
in the the magnetization curve m = m(h) at a criti- 
cal field h c w 0.941, m c = m(h c ) 0.556. In the low 
field regime h < h c , all three probabilities P + = N+/N, 
Pa = N /N and P_ = N-/N tend to finite values as 
N — * oo, but for h > h c the ground state is in a sector 
with = 0, i.e., P_ = 0. Thus P_ can be used as an 
order parameter for the transition. The magnetization 
behaves continuously at the critical point, but the sus- 
ceptibility diverges below h c . The high- field phase (SI 
phase), where the low-energy sector is identical to that 
of the spin-1/2 Heisenberg chain, can be described by a 
one-component Luttinger liquid, and there is only one 
characteristic soft mode. The low field phase (S2 phase), 
on the other hand, is equivalent to a two-component Lut- 
tinger liquid, possessing two critical degrees of freedom, 
and thus two incommensurate soft modes. Their posi- 
tions in the Brillouin zone are functions of the probabil- 
ities -P+,0.-- In the sequel, we label phases by Sn, where 
n = 0,1,2 stands for the number of critical, Luttinger 
liquid components. Note that SO is a gapped phase. 
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Early speculations™] that the S2-S1 phase transition 
of the ULS model may also take place in the Heisen- 
berg chain at f3 = was finally refuted by Takahashi and 
SakaiB who found that in the whole range < m < 1 the 
low-energy physics is described by a c = 1 U(l) confor- 
mal field theory (CFT) which is equivalent to the one- 
component Luttinger liquid, thus only an SI phase ap- 
pears. The Luttinger liquid parameters vary smoothly 
as a function of the magnetization. Naturally arises the 
question whether the appearence of the Haldane gap in 
the ground state only allows the SI behavior seen for 
the pure Heisenberg model, or there is a certain domain 
in parameter space where multi-component Luttinger liq- 
uids such as an S2 phase can occur above the SO Haldane 
phase. To clarify this question is the principal aim of this 
Letter. 

The first indication that the Haldane gap of the 
bilinear-biquadratic model may collapse into an S2 phase 
comes from the numerical observation that at h = 
the VBS point /3vbs = 1/3, where the ground state 
can be constructed exactly using nearest-neighbor va- 
lence bonds, is in fact a disorder point, beyond which 
short range fluctuations in the ground state become 
incommensurate. □ However, due to the finite correla- 
tion length, the peak at ir in the static structure factor 
only splits somewhat later at the so called Lifshitz point 
Puts ~ 0. 438.0 One can define a third special point /?Disp> 
which is a priori distinct from (but close to) the above 
two, where the emerging incommensurability make the 
position of the minimum gap in the energy-momentum 
dispersion relation move away from it. In the range 
/3disp < P < 1 the momentum pn associated with this 
gap minimum rapidly shifts from the antiferromagnetic 
value 7r to the ULS value 2ir/3. When the magnetic field 
reaches the value of the gap at pn, and the system starts 
to become magnetized, there are obviously two character- 
istic momenta in the system: one is 2-7rm, as suggested by, 
the generalization of the Lieb-Schultz-Mattis theorem,tll 
and the other is the finite difference of the two split gap 
minima 2pu- Of course, any linear combinations of these 
two may also show up in the correlation functions, which 
will be dominated by the most slowly decaying terms, 
or on shorter distances, by the ones having the largest 
amplitude. 

A more quantitative analysis is possible on the basis of 
the Bethe Ansatz (BA) solution of the ULS model, when 
the deviation term proportional to (1 — (3) in the Hamil- 
tonian is treated as a perturbation. The ULS model is 
solvable by the two-component nested BA method.B This 
associates the + spin states with an inert background in 
which particles with two possible internal states, spin 
and — , move. Two sets of spectral parameters are intro- 
duced, one for the particles (component-1) and one for 
their internal state (component-2). Their actual values 
can be calculated by solving a set of coupled algebraic 
(or in the N — > oo limit integral) equations. The BA 
technique also allows one to obtain the finite size cor- 
rections 0(1/N) to the low-energy excitations near the 
thermodynamic limit. In the S2 phase of the ULS model 



the energy spectrum has the following structured 
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where E g is the energy of the ground state, the in- 
dex a = {di, d2',h,l2'i tt-i, , , n^} is a shorthand for 
eight integer (half-integer) quantum numbers specifying 
the eigenstate, and the matrix Z a p, a, f3 = 1,2 is the 
"dressed charge matrix" responsible for the interaction 
of the two BA components. The relative momentum of 
the state a reads 



2tt 



SP = P a - P g = Q + jy(Ai - A^ + A^ - 
where Q = Q a is an 0(1) term 

Q = 2tt(1 - P+)di + 2irP-d 2 + irh. 



(4) 
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The physical interpretation of the quantum numbers is as 
follows: I a (d a ) represents the number of particles added 
to (transferred from the left Fermi point to the right in) 
component a. is the number of small momentum 
particle-hole pairs created in component a around the 
left (— ) and the right (+) Fermi points. While l a and 
n a — are always integers_d Q is integer or half integer 
with di t 2 = h,i/2 (mod This structure is analogous 
to the one found in the ID Hubbard model where the 
two components are called "charge" and "spin", resp. In 
the present case l\ = SNo + 5N-, I2 = 8N- for which 
there are selection rules when a given type of correlation 
functions is considered. 

The low-energy excitations of the ULS model in its 
S2 phase can be interpreted by assuming that they arise 
from the direct sum of two c = 1 conformal field theories 
(CFT-), each having a different sound velocity v\ and V2, 
resp.E3 As is indicated by Eq. (^) local physical operators 
necessarily couple to both CFTs. Conformal invafiance 
then requires that the 2-point functions behave agl3 



(<j>(x, 0)) = A^e~ lQx ]J(xt iv a t) 
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showing the analog of "spin-charge separation" for the 
present spin-1 situation. Let us consider the operator 
content of the theory: each operator </> a (primary or sec- 
ondary) is labeled by the eight quantum numbers a, and 
has the anomalous dimension x a — At + Ar + A^ + A, 
and conformal spin s a 
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A+ 



"1 ' ^2 ^2 — 
Note that the total mo- 



dih +d 2 l2 + n{ —1^ 
mentum associated to these operators is SP in Eq. (Q), 
involving the Q term as well. There are four marginal op- 
erators M. 1,2,3,4 which can be formed using only the n 1 2 
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quantum numbers and setting l a ,d a = 0. Their dimen- 
sion, spin and total momentum (x = 2, s = 0, SP = 0) 
do not depend on the Z a p matrix. The presence of these 
operators causes the existence of an extended critical 
domain in the parameter space. The other relevant or 
marginal operators are all primary, i.e., = and all 
depend on the Z a p matrix. This latter can be calcu- 
lated numerically for the ULS model by|-aolving a set of 
coupled integral equations. The resultsEJ are shown in 

Fig. |(a). When h = 0, Z n = Z 22 = yJl/3 + 1/2^3, 

Z12 = Z21 = yl/3~ — 1/2 V3, and there is an additional 

marginal operator </>i/2,i/2;i,-i (and its equivalents un- 
der SU(3) and parity transformations, e.g., </>i,i) as can 
be checked using Eq. (m. (From now on in the index 
a we omit l a and/or when they are zero.) The op- 
erator 0i/2,i/2 ; i,-i has anomalous dimension x = 2 and 
total momentum SP = when h = 0. This is the princi- 
pal operator responsible for the SU(3)^SU(2) symmetry 
breaking processes 00 < — ► + — , — h As was shown in 
Ref.cl the interplay of <^i/2,i/2;i,-i an d the M. operators 
(which constitute the SU(3) current interaction in the 
WZWN theory) gives rise to the Haldane gap for /3 < 1 
but maintains criticality for [3 > 1. 
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FIG. 1. Bethe Ansatz results for the ULS model vs the 
magnetization m. (a) The dressed charge matrix, (b) anoma- 
lous dimension and (c) momentum Q of some selected opera- 
tors. Some other operators marginal at h — are not shown, 
(d) DMRG results for the Fourier transform of the one-point 
function (S^). The peak at the Lieb-Schultz-Mattis value 
q = 2nm, 2ir(l — m) only develops for m > m c . 

When < h < h c the anomalous dimensions and mo- 
menta of the operators present in the ULS model vary as 
shown in Fig. [l](b). Since the perturbation part of the 
Hamiltonian, represented by the deviation (1 — 0) term, 
transforms under translations with momentum zero, in 
its decomposition into the operators present in the ULS 
model at h > only operators with SP = can appear. 



This is a serious limitation, since as shown in Fig. [l](c), 
the two characteristic momenta Qi,o = 27r(l — P + ) and 
Qo,i = 27rP_, associated to the large momentum trans- 
fer processes of the two components in Eq. (|^), become 
generically incommensurate. </>i/2,i/2;i,-i has no longer 
SP = so it does not contribute. It is in fact an Umk- 
lapp operator which appears in the low-energy descrip- 
tion only at h = 0. What contribute are the operators 
^di d 2 -i -i-n ± -n ± W1 th d\ , d,2 half- integer, and > cho- 
sen in a way to reinstall SP = 0. However, due to the 
appearance of the necessary small momentum particle- 
hole excitations such operators are highly irrelevant. The 
marginal operators M. are not able alone to drive the 
two Luttinger liquid components away from criticality, 
although they make the universality class change con- 
tinuously. We conclude that when < m < m c the 
bilinear-biquadratic model must remain in its S2 phase in 
an extended domain on both sides of the ULS line (3—1. 
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FIG. 2. Phase diagram of the bilinear-biquadratic spin-1 
chain in a magnetic field. The one- and two-component 
Luttinger liquid phases are denoted by SI and S2, resp. 
The S1-S2 phase boundary is determined by the DMRG (O 
points) . The dotted line indicates some uncertainties near the 
disorder points. 

When the magnetic field reaches the critical value h c 
in the ULS model the sound velocity V2 — » 0, and the 
corresponding critical degree of freedom becomes mas- 
sive. The emerging SI phase can be described by a single 
c = 1 CFT, and the universality class is determined by 
a scalar dressed charge Z. Once again we do not expect 
any drastical changes in the low-energy physics as (3 is 
perturbed away from 1. The predicted phase diagram is 
shown in Fig. |^. 

In order to find the exact phase boundaries we car- 
ried out a detailed numerical analysis using the DMRG 
technique. Unfortunately the DMRG does not work in 
momentum space thus a direct search for tracking soft 
modes by calculating energy gaps is not feasible. In- 
stead, we calculated the decay of the one-point function 
(Sj) from the edge of an open chain for different val- 
ues of (3 and magnetization m = S* ot /N. The one-point 
function contain the same information as the equal time 
two-point correlation function in Eq. (^J) (note that the 
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exponent of the one-point function is half of that of the 
two-point function), and the appearing soft modes can 
be identified in the Fourier transforms 0|t-.by making a 
suitable multi-parameter fit in real space.E3 

Once again exact results can be obtained for the ULS 
model. Considering (S?) in the S2 phase only operators 
with SS* ot — 0, i.e., h = 1% = and h = —h = ±1 can 
contribute. The most relevant operators are </>i.o, <Ao,i, 
0i/2,i/2;i,-i) and 0i/2,-i/2;i,-i (and their equiv- 
alents under symmetries). The associated momenta Q 
(position of the peak in the structure factor) and the oc- 
curring critical exponents can be read off from Fig. |l](b) 
and (c). The example presented in Fig. |l|(d) illustrates 
that in the S2 phase of the ULS model only the first three 
with l\ = 1% = have nonvanishing amplitude as dic- 
tated by the higher symmetry. However, 0i/2,±i/2;i,-i 
contributes when /3 ^ 1. It is remarkable that in the 
S2 phase the operator 0i,i, which has an anomalous di- 
mension over 2 and momentum Qi^i = 2tv(1 — to), does 
not contribute to the correlation function. In the high- 
field SI phase, however, the only peak discernible in the 
structure factor, as shown in Fig. |l](d), is the one with 
Qi = 2it(1 — rriL in agreement with the Lieb-Schultz- 
Mattis theorem O The presence or absence of a peak in 
the structure factor at k = ±2-7tto can thus be used effi- 
ciently to distinguish between the two phases and locate 
the phase boundary. Alternatively, one can monitor the 
amplitude of the peak at Qo,i which tends to zero on the 
transition line where the multi-peak structure collapses 
into a single-peak one. The phase diagram, as determined 
by the DMRG calculation, is shown in Fig. |^. Details pf, 
the numerical investigation will be published elsewhere.Ej 

Our analysis is valid in the strict sense close to the 
ULS line only. Here the S2-S1 phase transition is clearly 
associated to the depletion of one of the two "bands" 
as suggested by the BA. Note that in Fig. || the phase 
boundary does not change very much until about (3 ~ 

0. 5 where it seems to decline rather rapidly. A priori 
we cannot exclude the possibility that some operators 
become relevant here and open a gap in one of the critical 
components. This question needs further clarification. 

In the bilinear-biquadratic model in Eq. (|l|) the S2 
phase terminates near j3 » 0.4. This is still far in the 
parameter space from the currently known spin-1 Hal- 
dane gap materials for which the biquadratic term is 
small, and the pure Heisenberg model, although with 
some anisotropies, is a good description. For these sys- 
tems only SI type massless phases (and eventually, for 
some special values of to,, additional SO type phases, 

1. e., magnetization plateausElHU) are expected to appear. 
However, even here, due to the closeness of the S2 phase 
somewhat further in the phase diagram, massive but rel- 
atively low energy excitations are predicted to show up 
in the weakly magnetized regime. They are expected to 
contribute in experiments probing higher lying excita- 
tions such as in inelastic neutron scattering, or in situ- 
ations when short distance physics is important as, e.g., 
in nonmagnetically doped materials. Although there is 
no sharp phase transition in this case, the low-field and 



high- field regimes may look rather different, separated by 
a more or less narrow crossover region as observed, e.g., 
in Ref. 7. 

In general, a two-component Luttinger liquid phase 
(S2), and then an eventual phase transition S2— >S1 dur- 
ing the magnetization process, is expected to occur when- 
ever the ground state develops incommensurate fluctua- 
tions already at h = 0. It must not necessarily be above 
a Haldane gap; the spin- 1/2 zig-zag ladder, for exam- 
ple, which is expected to describe adequately the quasi- 
1D antiferromagnet CS2CUCI4, where the gap is due to 
dimerization and fluctuations are also predicted tp-.be 
incommensurate already without a magnetic fieldjij is 
another possible candidate. 
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